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ABSTRACT 

The quest for a 5-mode imprint from primordial gravity waves on the polarization of the cosmic microwave background (CMB) requires the 
characterization of foreground polarization from Galactic dust. We present a statistical study of the filamentary structure of the 353 GHz Planck 
Stokes maps at high Galactic latitude, relevant to the study of dust emission as a polarized foreground to the CMB. We filter the intensity and 
polarization maps to isolate filaments in the range of angular scales where the power asymmetry between E-modes and 5-modes is observed. 
Using the Smoothed Hessian Major Axis Filament Finder (SMAFF), we identify 259 filaments at high Galactic latitude, with lengths larger or equal 
to 2° (corresponding to 3.5pc in length for a typical distance of 100pc). These filaments show a preferred orientation parallel to the magnetic 
field projected onto the plane of the sky, derived from their polarization angles. We present mean maps of the filaments in Stokes /, Q , U, E, and 
5, computed by stacking individual images rotated to align the orientations of the filaments. Combining the stacked images and the histogram of 
relative orientations, we estimate the mean polarization fraction of the filaments to be 11 %. Furthermore, we show that the correlation between the 
filaments and the magnetic field orientations may account for the E and 5 asymmetry and the C EE /C EE ratio, reported in the power spectra analysis 
of the Planck 353 GHz polarization maps. Future models of the dust foreground for CMB polarization studies will need to take into account the 
observed correlation between the dust polarization and the structure of interstellar matter. 
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1. Introduction 

Recently, Planck 1 has reported an asymmetry in power between 
the dust E- and 5-modes in its 353 GHz observations (Planck 
Collaboration Int. XXX 2014; Planck Collaboration X 2015). 
This power asymmetry has been observed outside masks cov¬ 
ering 20 to 70 % of the sky, excluding the Galactic plane. The 
ratio of the dust 5- to 5-mode power amplitudes is about a half 
over the multipole range 40 < i < 600 (Planck Collaboration 
Int. XXX 2014). The source of this power asymmetry in the dust 
polarization data is currently unknown. Models of the Galactic 
magnetic field (GMF) used in the Planck sky model (PSM, 
Delabrouille et al. 2013) and the FGPol model (O’Dea et al. 
2012) produce an equal amount of power in E- and 5-modes 
outside the regions covered by the sky masks. These models, 
which were used to estimate the dust polarization foreground 
(BICEP2 Collaboration 2014; BICEP2/Keck Array and Planck 
Collaborations 2015), include an analytical model of the large- 
scale GMF (5 G mf) anc * a statistical description of the turbulent 
component of the magnetic field. 

The Planck maps of thermal dust emission display filaments 
distributed over the whole sky (Planck Collaboration XI 2014). 
The filamentary structure of the diffuse interstellar matter is also 
a striking feature of dust observations at higher angular reso¬ 
lution, performed by Herschel , and of spectroscopic Hi obser¬ 
vations (e.g. Miville-Deschenes et al. 2010; Andre et al. 2014; 
Clark et al. 2014). The analysis of Planck dust polarization data 
in the diffuse interstellar medium (ISM), at low and intermedi¬ 
ate Galactic latitudes, indicates that the structures of interstellar 
matter tend to be aligned with the plane of the sky (POS) pro¬ 
jection of the magnetic field (2?pos> Planck Collaboration Int. 
XXXII 2014). This preferential relative orientation is also ob¬ 
served in simulations of magneto-hydrodynamic (MHD) turbu¬ 
lence of the diffuse ISM (Hennebelle 2013; Soler et al. 2013). 
Such a coupling between the structure of interstellar matter and 
2?pos is not included in the PSM or FGPol models of the dust 
polarization sky (Delabrouille et al. 2013; O’Dea et al. 2012). 

The goal of this paper is to test whether the correlation be¬ 
tween the filamentary structures of the intensity map and 5 P os 
in the diffuse ISM accounts for the observed E-B asymmetry. 
Zaldarriaga (2001) describes the E- and 5-modes decomposi¬ 
tion of simple patterns of polarized emission, including filaments 
with a homogeneous polarization degree and orientation. The 
presence of 5-modes is related to invariance by parity of the po¬ 
larization pattern. There is 5-only power if 5pos is either paral¬ 
lel or perpendicular to the filaments. If 5 P os is oriented at +45° 
or -45° with respect to the filaments, there is 5-only power. 

In this paper, we filter the Planck intensity and polarization 
maps to isolate filaments in the range of angular scales where the 
5-5 asymmetry is observed. We identify coherent elongated fila¬ 
ments within regions of low column density at high Galactic lat¬ 
itude using a filament-finding algorithm. We evaluate the mean 
polarization angle in each of these filaments and compare it to 
the mean orientation of each filament. In doing so, we extend the 
analysis presented in Planck Collaboration Int. XXXII (2014) to 
the relevant region of the sky for CMB polarization observations 
at high Galactic latitude. In order to achieve a high signal-to- 


1 Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci¬ 
entific consortia funded by ESA member states and led by Principal 
Investigators from France and Italy, telescope reflectors provided 
through a collaboration between ESA and a scientific consortium led 
and funded by Denmark, and additional contributions from NASA 
(USA). 


noise ratio and enhance the contrast with respect to the local 
background dust emission, we stack the Stokes /, Q , U , and also 
5 and 5 maps, for the filaments we select in the Planck dust in¬ 
tensity map. We use the stacked images to quantify the power 
asymmetry in 5- and 5-modes associated with the filaments. 

This paper is organized as follows. In Sect. 2, we introduce 
the Planck 353 GHz data used in this study. The filament-finding 
algorithm is presented in Sect. 3. Section 4 presents the study of 
relative orientation between the filaments, 5 P0 s, and the POS 
component of the large-scale GMF (5 m? pos) at high Galactic 
latitude. In Sect. 5, we present the stacking of both intensity 
and polarization maps and derive the mean polarization frac¬ 
tion of the filaments. In Sect. 6, we discuss the relation between 
the relative orientation of the filaments and 5 P0 s and the 5-5 
asymmetry. Section 7 presents our results in the context of ear¬ 
lier studies and its relation to Galactic astrophysics. Finally, we 
present our conclusions in Sect. 8. This paper has four appen¬ 
dices. Appendix A details the Hessian analysis implemented to 
identify the filaments in the dust intensity map. The application 
of the filament-finding algorithm to a simulated Gaussian dust 
sky is detailed in Appendix B. In Appendix C, we study the im¬ 
pact of our selection of the filaments on the main results of the 
paper. The computation of all the angle uncertainties that we use 
in our analysis is presented in Appendix D. 


2. Planck data 

The Planck satellite has observed the sky polarization in seven 
frequency bands from 30 to 353 GHz (Planck Collaboration I 
2015). In this paper, we only use the 2015 (“DXlld”) data from 
the High Frequency Instrument (HFI, Lamarre et al. 2010) at 
353 GHz, since they are best suited to study the structure of 
dust polarization (Planck Collaboration Int. XIX 2015; Planck 
Collaboration Int. XX 2015; Planck Collaboration Int. XXI 
2015; Planck Collaboration Int. XXII 2015). The data process¬ 
ing, map-making, and calibration of the HFI data are described 
in Planck Collaboration VII (2015) and Planck Collaboration 
VIII (2015). In our analysis, we ignore the dust and CO spec¬ 
tral mismatch leakage from intensity to polarization (Planck 
Collaboration VIII 2015). Planck Collaboration Int. XXX (2014) 
has shown that the amplitude of the dust spectral mismatch leak¬ 
age at high latitude (/ s k y = 0.5) is small compared to the total 
polarization signal in 5 and 5 modes. No CO emission is de¬ 
tected at 353 GHz away from the Galactic plane and the bright¬ 
est molecular clouds (Planck Collaboration XIII 2014) and so 
we do not consider it in our analysis. 

To quantify the statistical noise and systematic effects on 
the results presented in this paper, we use the two HalfMission 
(HM), two DetSet (DS), and two HalfRing (HR) Planck 
353 GHz polarization maps (Planck Collaboration VII 2015). 
The two HM maps are made from the two halves of the full- 
mission Planck data, whereas the two HR maps are produced by 
splitting each ring (also called stable pointing period) into two 
equal duration parts. The two DS maps are constructed using 
two subsets of polarization-sensitive bolometers at a given fre¬ 
quency. The noise is uncorrelated between the two HM, HR, and 
DS maps. We only use them to compute the error bars on the 
relevant quantities that we measure in this paper. 

The total polarization intensity (P 353 ) and the polarization 
angle (if/) are derived from the full-mission Stokes Q 353 and U 353 
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maps at 353 GHz using the relations 

^353 = ^6353 + £353 ? ( 1 ) 

if/ = 0.5 x atan 2 (-f/ 353 , £> 353 ), (2) 

where the two-argument function atan 2 (-£ 353 , ^ 353 ) is used to 
compute atan(-[ 7353 / 2353 ) avoiding the n ambiguity. To recover 
the correct full range of polarization angles ([—tt/2, 7t/2 ] as used 
for if/ here), attention must be paid to the signs of both £353 
and 2353? not just their ratio. We use the IAU convention for 
if'f, which is measured from the Galactic North (GN) and posi¬ 
tive to the East. The minus sign in Eq. (2) converts the conven¬ 
tion provided in the Planck data to that of the IAU (see Planck 
Collaboration Int. XIX 2015). The orientation angle (X) of £ P0 s 
is defined within the n ambiguity by adding zr /2 to the polariza¬ 
tion angle 

X = + 2 • (3) 

For the dust intensity at 353 GHz, we use the model map 
7 ) 353 , computed from a modified blackbody fit to the Planck 
data at v > 353 GHz and the IRAS 100 jam map (Planck 
Collaboration XI 2014). This map has lower noise than the 
353 GHz Stokes 7 map and is corrected for zodiacal light emis¬ 
sion, CMB anisotropies, and the cosmic infrared background 
monopole. We neglect the contribution of the CMB polarization 
at 353 GHz for this study. 

The full-mission Planck Stokes 2353 and £353 maps are pro¬ 
vided in HEALPix 2 format (Gorski et al. 2005) at 4 .'8 resolu¬ 
tion and D 353 at 5'. To increase the signal-to-noise ratio, we 
smooth the three maps to a common resolution of 15', taking 
into account the effective beam response of each map, and re¬ 
duce to a HEALPix resolution of N s id e = 512. For the polariza¬ 
tion data, we decompose the Stokes 2353 and £353 maps into 
£■353 and £353 ai m s (£^ m and £^ m ) using the “ianafast” routine 
of HEALPix, apply the Gaussian smoothing in harmonic space 
(after deconvolving the effective azimuthally symmetric beam 
response of each map), and transform the smoothed £353 and 
#353 atm s back to 2353 and £353 maps using the “isynfast” 
routine at = 512. We also transform the £353 and £353 at m s 
to £353 and £353 maps at = 512 using the relations 

£353(11) = ^ EimYimin ), £353(11) = ^ B{ m Yi m (n) . (4) 

All the maps that we use are in thermodynamic units OuKcmb)- 

In this paper, we work with the bandpass-filtered dust in¬ 
tensity map, 7)!^, to identify and isolate filaments over the fil¬ 
tering scale using a filament-finding algorithm. By filtering out 
large-scale and small-scale modes, we enhance the contrast of 
the filaments with respect to the diffuse background and reduce 
the instrumental noise, which is critical for accurately measur¬ 
ing the polarization orientations of the filaments within regions 
of low column density at high Galactic latitude. 

For filtering, we apply the three-dimensional spline wavelet 
decomposition based on the undecimated wavelet transform, as 
described by Starck et al. (2006). We use the publicly avail¬ 
able package Interactive Sparse Astronomical Data Analysis 
Packages (ISAP 3 ) to compute the D ^ 53 map at A S id e = 512 reso¬ 
lution. The spline wavelet used in this analysis provides less os¬ 
cillation in position space compared to Meyers or needlet ones 
(Fanusse et al. 2012). The filtering is done in pixel space; the 
corresponding bandpass filter in harmonic space varies a little 

2 http://healpix.jpl.nasa.gov 

3 http://www.cosmostat.org/isap.html 
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Fig. 1 . The representative bandpass filter, retaining only the 
scales between l - 30 and 300. 


over the sky. Figure 1 presents the typical shape of this band¬ 
pass filter, which selects the scales between i - 30 and 300. The 
filtering scale is chosen in such a way that it highlights all the 
bright filaments present in the Planck 7)353 map. We also com¬ 
pute the bandpass-filtered polarization maps, 2353 ? £ 353 ? £ 353 ’ 
and £ 353 . 

3. Filament-finding algorithm 

3.1. Methodology 

Identification of filaments as coherent structures is a cru¬ 
cial part of this analysis. Previous studies used algorithms 
such as DisPerSE (Sousbie 2011; Arzoumanian et al. 2011), 
getfilaments (Men’shchikov 2013), and the rolling Hough 
transform (Clark et al. 2014). Hennebelle (2013) and Soler et al. 
(2013) used the inertia matrix and the gradient of the density and 
column density fields to identify filaments in numerical simula¬ 
tions of MHD turbulence. 

In this paper, we employ the Smoothed Hessian Major Axis 
Filament Finder (SMAFF, Bond et al. 2010a) algorithm, which 
has been used to identify filaments in the three-dimensional 
galaxy distribution (Bond et al. 2010b). SMAFF is primarily based 
on the Hessian analysis. The Hessian analysis has also been 
used to analyse the Planck dust total intensity map in Planck 
Collaboration Int. XXXII (2014), Herschel images of the FI641 
cloud in Orion A (Polychroni et al. 2013), and large-scale struc¬ 
ture in simulations of the cosmic web (Colombi et al. 2000; 
Forero-Romero et al. 2009). Planck Collaboration Int. XXXII 
(2014) has reported good agreement between the filament orien¬ 
tations derived from the Hessian and inertia matrix algorithms. 

3.2. Implementation 

In this study, we apply the two-dimensional version of SMAFF 
to the 7)353 map, which is shown in the upper right panel of 
Fig. 2. From the Hessian matrix, we compute an all-sky map 
of the lower eigenvalue /L and the orientation angle 6 - of the 
perpendicular to the corresponding eigenvector, measured with 
respect to the GN. The details of the Hessian analysis are pro¬ 
vided in Appendix A. The map of /L is presented in the lower 
left panel of Fig. 2. For the subsequent analysis, we consider 
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Fig. 2. Data processing steps implemented to identify filaments from the Planck data. We start with the Planck D 353 map (upper 
left panel) smoothed at 15' resolution. The bandpass-filtered D 353 map (upper right panel) is produced using the spline wavelet 
decomposition, retaining only the scales between € = 30 and 300. The lower eigenvalue map of the Hessian matrix, T_, is shown 
in the lower left panel. Structures identified in the high-latitude sky A- map are shown in the lower right panel. The superimposed 
graticule is plotted in each image and labelled only on the lower right panel. It shows lines of constant longitude separated by 60° 
and lines of constant latitude separated by 30°. The same graticule is used in all plots of the paper. 


only the high-latitude sky, defined as \b\ > 30°, with the Large 
Magellanic Cloud and Small Magellanic Cloud regions masked 
out. 

The map of A- highlights filaments in the D 353 map with 
an orientation angle 6_, which we refer to as the Hessian an¬ 
gle hereafter. The distribution of /L over the unmasked pixels 
is shown in Fig. 3. This distribution of A- is non-Gaussian with 
an extended tail. We use the median absolute deviation (MAD, 
Hampel 1974; Komm et al. 1999) to measure the width, <x /L , of 
the distribution, as given by 

crx_ - 1.4826 x median(|T_ - m^J), (5) 

where m / t_ is the median of the A- distribution. We select fila¬ 
ments using an upper threshold (K_) on T_ given by 

K- = mx_ - 3 cr A _ . ( 6 ) 

Hereafter, we refer to the filaments satisfying A- < K- as 
“strong”. This threshold K- separates the strong filaments from 
the weak ones, as detailed in Appendices B and C. By construc¬ 
tion, the threshold K- rejects pixels where A- is positive, since 
those pixels do not correspond to local maxima. 

We seek coherent elongated structures in the map. In SMAFF, 
this is achieved by placing an upper limit C on the difference 



A- [K C mb rad 2 ] 

Fig. 3. Distribution of the eigenvalues /L over the unmasked pix¬ 
els in the high-latitude sky. The grey region represents the pixels 
that were used in the SMAFF algorithm to find strong filaments. 
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between Hessian angles within a given structure. For our pur¬ 
pose, we set the value of C = 15° to identify relatively straight 
filaments. 

We start with the pixel having the most negative A- and 
denote the corresponding Hessian orientation angle by 0 S _- We 
identify its neighbouring pixels using the “neighbours_ring” 
routine of HEALPix and look for pixels with A- < K- and orien¬ 
tation angle such that 10_ - 6 S _ | < C. If both conditions are sat¬ 
isfied, we count that neighbouring pixel as a part of the filament 
and move on to that neighbouring pixel. The neighbouring pixel 
becomes the new reference point and we search for its neigh¬ 
bours that satisfy both conditions. In our algorithm, 6 S _ is fixed 
by the starting pixel, which has the most negative T_. We con¬ 
tinue this friend-of-friend algorithm to connect pixels until one 
of the conditions is no longer satisfied. We limit our selection 
to filaments with a length (L, defined as the maximum angu¬ 
lar distance between pixels within a given structure) larger than 
or equal to the threshold length Lo, which we choose to be 2 °. 
This process yields a set of 259 elongated filaments, as shown 
in the lower left panel of Fig. 2. Hereafter, we refer to this set as 
our filament sample. Selected sky pixels represent 2.2 % of the 
high-latitude sky considered in our analysis. There is no overlap 
between the filaments in our sample. 

The column density is computed from the D 353 map using 
the conversion factor, 0.039 MJysr -1 per 10 20 Hcm -2 ; this was 
derived in Planck Collaboration Int. XVII (2014) by correlating 
the Planck 353 GHz dust total emission map with an Hi column 
density map over the southern Galactic polar cap. We average 
the column density along each filament and assign one mean 
column density, A/^, to each. This column density is computed 
on the filtered intensity map. The histogram of for the fila¬ 
ment sample is presented in Fig. 4. The number of filaments per 
N ^ is represented by V F . 



A/b [10 20 cm -2 ] 

Fig. 4. Histogram of the mean column density of the filament 
sample. The column density is computed from the D 353 map us¬ 
ing the conversion factor derived in Planck Collaboration Int. 
XVII (2014). 



Fig. 5. Sketch of the mean orientation angle of the filament (6_), 
the magnetic field (?), the polarization angle (if/), and the large- 
scale GMF (X m ) along the filament. All the angles are defined 
with respect to the GN and follow the IAU convention. 


4. Interplay between the filament orientation and 
the magnetic field 

In this section, we study the orientations of matter structures and 
Bpos in our filament sample (Sect. 3). The orientation angle of 
Z?pos is derived from the observed Stokes £>353 and t /353 maps 
using Eqs. (2) and (3). We also consider the orientation angle 
(X m ) of # m ,POS> as estimated from starlight polarization observa¬ 
tions (Heiles 1996) and pulsar rotation measures (Rand & Lyne 
1994; Han et al. 1999). We compare these three orientations, as 
represented in Fig. 5. Our analysis follows Planck Collaboration 
Int. XXXII (2014), which used a set of pixels representing ap¬ 
proximately 4 % of the sky at low and intermediate Galactic lat¬ 
itudes. Only 25 % of the pixels in our current filament sample 
were considered in this earlier study. 

4.1. Relative orientation of the filaments and the magnetic 
field 

We study the angle difference between the orientations of the fil¬ 
aments and Bpos in our sample. First we associate one POS ori¬ 
entation angle with each of the filaments with respect to the GN. 
By construction, due to our selection criteria on the angles, the 
filaments are fairly straight and, hence, they may be described 
with a single orientation angle. Given one filament, we measure 
the mean orientation angle, 0 _, over the n pixels that belongs to 
it. We make use of the pseudo-vector field with unit length com¬ 
puted from the values of 0_ for each pixel. This pseudo-vector 
has components g- = cos 20 _ and U- = - sin 26- (following 
the HEALPix convention for the g_ and U- components). The 
mean POS orientation angle 0_ of the filament is obtained by 
first averaging g_ and U- over all n pixels and then calculating 
the position angle of this averaged pseudo-vector. It is given by 

( 1 n , n \ 

— Y U-,-Y Q- . (7) 

n U ) 
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If we rotate the Stokes £>353 and C /353 maps by 0_, i.e., into the 
frame where the axis of the filament is in the North-South direc¬ 
tion, the rotated 2353 and C /353 can be written as 

6353 = 2353 cos 26 _ - U 353 sin 26- , ( 8 ) 

U 333 = 2353 sin 26- + U 333 cos 26- . (9) 

Combining Eqs. ( 8 ) and (9) with Eqs. (1) and (2), we get 

2353 - ^353 cos 2 (if/ - 0_) = - P 353 cos 2(X - 6 L), (10) 

^353 = - P 353 sin 2 (iff - 6 -) = P 353 sin 2(X - 6 -) , (11) 

where the orientation angled is defined in Eq. (3). Similar to the 
computation of 6 -, we average Q' 353 and U' 353 over all n pixels 
and then calculate the position angle of this averaged pseudo¬ 
vector 

1 n 

2353 = — 2353 = — ^353 cos 2 Ax_e_ , ( 12 ) 

n i =1 

1 n 

U 333 = — ^353 = ^353 sin 2 A x-e_ » (13) 

n i =1 

where 

Ar- 0 _ = 0.5 x atan 2 (f? 353 , “ 2353 ) > (14) 

^353 = ^2 3 53 + ^353 • (15) 

The angle difference measures the weighted mean of the 
angle difference per pixel between the orientations of the given 
filament and Z? P 0 s- The index X refers to the mean orientation 
angle of Z? P0 s along the filament. Note that we directly mea¬ 
sure the angle difference between the filament and 5 P os, without 
computing X for each filament. 

The histogram of relative orientation (HRO) between the fil¬ 
ament and 2 ? P os for our filament sample is presented in the upper 
panel of Fig. 6 . The mean value of the histogram is 2?3 com¬ 
puted using the equivalent of Eq. (7). Our histogram agrees with 
the pixel-by-pixel analysis at intermediate and low Galactic lati¬ 
tudes presented in Planck Collaboration Int. XXXII (2014). Like 
in this earlier study, we find that the filaments are statistically 
aligned with Bp os- A similar alignment between the filaments 
in the intensity map and 2 ? P os has been reported for synchrotron 
emission observed by WMAP at 23 GHz (Vidal et al. 2015). 

To quantify the shape of the histogram of A x-e _, we fit it with 
a Gaussian plus a constant. The Gaussian has a 1 cr dispersion of 
19°. The constant may be accounted for by the projection of the 
magnetic field and filament orientations on the POS as discussed 
in Planck Collaboration Int. XXXII (2014). 

4 . 2 . Relative orientation of the magnetic field and the 
large-scale Galactic magnetic field 

Here, we compare the orientation of 2? P os on the filaments with 
that of # m ,Pos- Heiles (1996) derived the orientation of #m,Pos 
pointing towards Z 0 = 82?8 ± 4?1 and bo = 0?4 ± 0?5 from 
the polarization pseudo-vectors of stars more distant than 500 
parsecs. Slightly different Zq values have been reported in other 
studies. From the rotation measures of nearby pulsars within a 
few hundred parsecs of the Sun, Rand & Lyne (1994) found the 
direction of # m ,POS pointing towards Zo = 88 ° ±5°. In another 
study of pulsar rotation measures, Han et al. (1999) derived the 
direction of B m pos as Zo - 82°. These two studies do not report 





Fig. 6. Upper panel : HRO between the filaments and 5 P os- 
Middle panel : HRO between 2Z P os an d # m ,POS- Lower panel : 
HRO between the filaments and # m ,Pos • 


values for Z?o, which is assumed to be zero. Based on these ob¬ 
servations, we assume that the mean orientation of # m , P os in the 
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Fig. 7. Mean images of the Planck (£>353) , (G353X and (t/' 53 > maps over the T_ filaments at 15' resolution. 



solar neighbourhood is Iq = 84° ±10° and bo = 0° ± 10°, with 
the same uncertainty on / 0 and bo. 

We construct a pseudo-vector field with unit length based 
on the uniform orientation of # m ,POS- This pseudo-vector has 
components: Q m = cos 2if/ m = cos 2(X m - n/2) and U m = 
- sin 2 if/ m = - sin 2(X m - 7 r/ 2 ) (following the HEALPix con¬ 
vention for the Q m and U m maps), where if/ m is the polarization 
angle of # m ,POS- The procedure to go from the uniform B m PO s 
pointing towards (/o, bo) to i// m is detailed by Heiles (1996). The 
mean orientation angle (X m ) of # m ,POS for each filament is ob¬ 
tained by first averaging Q m and U m over all n pixels within a 
filament and then calculating the position angle of this averaged 
pseudo-vector. We compute the angle difference, A^_^ m , between 
the orientations of 2 ?pos and Z? m ,POS on the filament in a similar 
manner to the method described in Sect. 4.1. 

The HRO between Z?pos and Z? m ,Pos for our filament sample 
is presented in the middle panel of Fig. 6 . The mean value of 
the histogram is 1?3±3?7 where the uncertainty is computed by 
changing the mean orientation of # m ,POS within its quoted un¬ 
certainties. This HRO has a larger dispersion than that between 
the filaments and 2 ?pos shown in the upper panel of Fig. 6 . To 
quantify the shape of the histogram of A^_^ m , we fit it with a 
Gaussian plus a constant. The Gaussian has a 1 cr dispersion of 
36°. 

We conclude that Z? P os of our filament sample is statistically 
aligned with 2? m ,POS- Planck Collaboration Int. XXXII (2014) re¬ 
ports a similar correlation, for the low and intermediate Galactic 
latitudes, when comparing the polarization measured on the fila¬ 
ments with their background polarization maps. The scatter mea¬ 
sured by the HRO may be interpreted considering both the tur¬ 
bulent component of the magnetic field and projection effects. 

4.3. Relative orientation of the large-scale Galactic magnetic 
field and the filaments 

We combine the results obtained in Sects. 4.1 and 4.2 to assess 
statistically the orientation of B m PO s in the solar neighbourhood 
with respect to the filaments. 2 ?pos is statistically aligned with 
the filaments in our sample and with # m ,POS- From both results, 
one would intuitively expect # m ,POS to be statistically aligned 
with the filaments. To test this expectation, for each filament, we 
compute the angle difference, A^ m _^_, between the orientations 
of # m ,POS and the filament. The angle difference A x m -o_ i s com¬ 
puted in a similar manner to the method described in Sect. 4.1. 


The HRO between the filament and # m ,POS for our filament 
sample is presented in the lower panel of Fig. 6 . A correlation 
between the orientation angles of the filament and Z? m ,POS is 
present, but the HRO shows more scatter than the HRO between 
the filaments and #pos and that between 5 P0 s and # m ,Pos- The 
histogram has a mean value of -3? 1 ± 2? 6 . To quantify the shape 
of the histogram of A^ m , we fit with a Gaussian plus a constant. 
The Gaussian has a 1 cr dispersion of 54°. Planck Collaboration 
Int. XXXII (2014) reported a similar loss of correlation when 
comparing the orientations of the filaments with that of 2 ?pos 
derived from their local background polarization maps. 

5. Mean polarization properties of filaments 

In this section, we present stacked images of the filaments in 
Stokes /, g, and U, after rotation to align the filaments and to 
compute Q and U with respect to their orientation. The images 
are used to compute the average polarization fraction of our sam¬ 
ple of filaments. 

5.1. Stacking filaments 

Over the high-latitude sky, the signal-to-noise ratio of the 
353 GHz Planck polarization maps is low and it is not possible 
to measure the polarization fraction of individual dust intensity 
filaments in our sample. In order to increase the signal-to-noise 
ratio, we therefore stack images of the 259 filaments and their 
surroundings. 

For each filament in the sample, using the “gnomview” rou¬ 
tine of HEALPix, we extract from the Planck maps a local, flat- 
sky, image (7° x 5° patch) centred on the filament centre and 
rotated by 6 - in the clockwise direction to align the filament in 
the North-South direction. We stack the images of the filaments 
in 7 ) 353 , 6353 ’ an d U' 353 (as defined in Eqs. 8 and 9) after align¬ 
ing all the maps in the North-South direction. We produce mean 
stacked images, denoted with angle brackets (..), by dividing the 
sum of the individual images by the total number of filaments in 
our sample; they are presented in Fig. 7. The 1 cr errorbar both on 
the ( Q 353 ) and (U' 353 ) images is 1.3 juK C mb, as computed from 
the difference of two polarization HM maps. All the features pre¬ 
sented in Fig. 7 are significant compared to the data systematics 
and statistical noise. The average filament appears as a negative 
feature with respect to the background in the (Q 353 ) image and is 
not seen in the (U ' 353 ) image. This result is a direct consequence 
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of the alignment between the filaments and Z? P os (Sect. 4.1 and 
Zaldarriaga 2001). The background in both the (Q 353 ) and ( U' 353 ) 
images is rather homogeneous. This reflects the smoothness of 
Z?pos within the 7° x 5° patches. 

We perform a null test to assess the significance of the stack¬ 
ing of filaments. This test is made by stacking 259 randomly 
chosen 7° x 5° patches in the high-latitude sky. Each patch is 
rotated in the clockwise direction, with the orientation angle 6 - 
of the central pixel. The images of (£> 353 ), (g 353 ), and (U ' 353 ) 
for random patches are consistent with noise. The amplitude of 
(Q 353 ) and (C/ 353 ) images is comparable to that of the difference 
between stacked images obtained when applying the same anal¬ 
ysis to each of the two polarization HM maps. This confirms the 
hypothesis that the filaments detected in Fig. 7 are indeed real 
and are rotated with a well-determined angle . 


5.2. Polarization fraction 


Instead of using individual pixels, we collapse the mean stacked 
images in the filament direction to draw the radial profiles ( R ) of 
the £>353, Q ' 353 , and U ' 353 images, which are presented in Fig. 8 . 
The shaded area in Fig. 8 represents the 1 <x dispersion from the 
data values at a given radial distance from the filament axis that 
we average. We clearly identify the profile of the filament on top 
of the constant background emission in the £>353 and Q 353 ra¬ 
dial profiles, while the radial profile of U ' 353 is consistent with a 
constant background emission. The radial profiles of £>353, Q 353 , 
and U 353 can be decomposed into the filament (F) and the back¬ 
ground (B) contributions as 

R d 353 = R d f + R d b , (16) 

333 ^353 ^353 

Rq' = Rn' F + Ro'Z 5 (17) 

^353 ^353 ^353 V ' 

Rjj' = Rtj'F + Rtj '-B • (18) 

U 353 U 353 U 353 V 7 

We fit the radial profiles of £>353 and Q 353 in Fig. 8 with a 
Gaussian profile for the filament emission plus a constant for the 
background emission. We find that the centre of the Gaussian 
profile is zero and that their 1 <x dispersion is 27' ± V for both 
the 7)353 and Q 353 radial profiles. 

Following Eqs. (12) and (13), we can express the average 
Stokes Q 353 and U ' 353 for one given filament as 

2353 = - /3 353 C0S 2A |_g_ = - f ^353 C0S 2A |-g_ > ( I9 ) 

^353 = ^353 S ' n 2 A J- 9 _ = / ^353 s ' n 2 A |Lg_ > ( 2 ®) 


where £) F 53 and P 353 are the average specific intensity and po¬ 
larization intensity of the filament. The superscript F represents 
the contribution from the filament only. The polarization fraction 
(p F ) of a filament is defined by 


pF 

r 353 

£> f 

^353 


( 21 ) 


From Fig. 4, we know that most filaments in our sample have 
comparable column densities and hence roughly the same £) F 53 . 
The mean stacked Stokes £2' F 3 and U ' 353 values for all the fila¬ 
ments can be approximated as 

< 6353 ) - -(/cos 2A|_ § )0 F 353 ) - -</><cos 2A|_ § )0 F 353 ), 

( 22 ) 

<^53) ^ </ sin 2A te >< 5 353> ^ </>< sin 2A |-g >< 5 353>’ 

(23) 
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Fig. 8. Radial profiles of the mean stacked Planck £> 353 , Q' 35V 
and U 353 images as functions of distance from the centre of the 
filament (blue line). The grey shaded region shows the 1 <x dis¬ 
persion from the data values at a given radial distance from the 
filament axis that we average. The dashed line is the Gaussian fit 
to the filament profile plus a constant background emission. 


where (p F ) is the mean polarization fraction of our filament sam¬ 
ple. For the radial profiles of the filament emission, we have 

Rq*, - -</><cos 2 A F _- g _)R D ? m , (24) 

Ru * 3 ^ </><sin 2 A |_g _>^ 53 • (25) 

The angle difference between the filaments and 2 ?pos (Eq. 14) 
is used for A F _- . The histogram of A F _- is roughly symmetric 
around 0°, implying Ru ,F 53 ^ Rq* 53 - 

Similarly, the radial profiles from the background emission 
can be written as 


R 

R 


o' B 

^ 35 : 


- -<p B ><cos 2 A|_ § _) j R d b 3 , (26) 

- (p B )(sin 2 A b _ 0 R d1:i , (27) 
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Rd 353 [mKcmb] 

Fig. 9. Linear fit (dashed line) to the correlation between the 
mean radial profiles of the stacked Q ' 353 and D 353 images. The 
grey shaded region shows the 1 cr dispersion from the data values 
at a given radial distance from the filament axis that we average. 

where (p B ) is the average polarization fraction of the back¬ 
ground emission and A?_- is the angle difference between a 
given filament and Z? P0 s from its local background polarization. 

The observed radial profile of the stacked Q 353 image in 
Fig. 8 can be written as 

Rq' = R Q ,F + Rrf B 

^353 ^353 Q 353 

- -<P F )<cos 2A|_ g _)7? D F 53 - (p B )(cos 2A|_ g _)/? D B 3 
= -<P F ><cos 2A?_- e _)R D353 

+ [</)<cos 2A|_ g _) - (p B )(cos 2A B _- e )\R Dl} 

— ci Rd 353 + b , (28) 

where a = -(/? F )(cos 2 A? - ) is the scaling parameter and 

b = [(p F )(cos 2 A|_^ ) - (p B )(cos 2 A?_^ )j R d b ^ is the offset of 
the linear fit between the Rq 353 and R D353 profiles. The linear fit 
between Rq 353 and R D353 is shown in Fig. 9. The best-fit parame¬ 
ter values from the linear fit are 

</><cos 2A|_g_) = 5.3%, (29) 

[</><cos 2 A|_ g ) - (p B >(cos 2 A|_ g )]r d ^ = 19.8/iKcmb . 

(30) 

In our estimate of (p F )(cos 2A|_- ), the noise bias is negli¬ 
gible, because the noise averages out in the stacking of the Q 353 
and U 353 maps. The HRO between the filament and 2 ?pos (upper 
panel in Fig. 6) is used to compute (cos 2A|_- ), which is 0.48. 
Taking this factor into account, the mean polarization fraction of 
our filament sample is 

( P F ) = n%. (3i) 

We have computed (p F ) using the two independent subsets of 
the Planck data (HM maps). For both data sets, we find the 
same mean value of 11 %. This is expected, because the HRO 
between the filaments and Z? P os is determined with high accu¬ 
racy (Appendix D). It shows that the measurement error on (p F ) 


is small. Obviously (p F ) may be different for another set of fil¬ 
aments, because it depends on the angles of the filaments with 
respect to the plane of the sky. 

The mean polarization fraction of the filaments determined 
from the filament sample is smaller than the maximum degree of 
polarization reported in Planck Collaboration Int. XIX (2015), 
which is /? max = 19.6 %. This result suggests that there is some 
depolarization due to changes in the magnetic field orientation 
within the filaments (Planck Collaboration Int. XX 2015; Planck 
Collaboration Int. XXXIII 2014). 

The offset b from the linear fit to the radial profiles of the 
Q 353 and D 353 images is positive. This means that 

<P B Xcos 2Af_g_) < </><cos 2A|_ g _) (32) 

The distribution of A B _- for our filament sample is needed to 

compute (p B ). We make a simple approximation using the HRO 
between the filament and B m pos (lower panel in Fig. 6 ) as a 
proxy for A ?_- Q . We then find (cos2A ?_- Q ) = 0.25 ± 0.14 for 
our filament sample. Combining Eqs. (29) and (32), we put a 
upper limit on the mean polarization fraction of the background 
emission as 

(p B )< 21%+?*. (33) 

This upper limit on (p B ) depends on the mean orientation of 
Z?gmf and it ranges from 48 % to 14 % for the 10° uncertainty 
on L and bn (Sect. 4.2). The profile of R r f B in the lower panel 

U 353 

of Fig. 8 is roughly constant and positive. This result indicates 
that Bpos is smooth within the 7° x5° patches and (sin 2A?_- ) is 

positive (Eq. 27), which follows if the distribution of A? - is not 
symmetric with respect to 0°. We point out that the histogram in 
the lower panel in Fig. 6 may not be used to estimate the sign of 
(sin2A?_- ) because the uncertainty on /o and bo is too large. 

6. E-B asymmetry 

In this section, we quantify the E-B asymmetry of dust polar¬ 
ization, over the l range 30 to 300, for our filament sample, and 
relate it to the relative orientation between filaments and 5 P os- 
First, we present the stacked images of filtered maps that we use 
to quantify the E-B asymmetry. Second, we compute the contri¬ 
bution of the pixels used in the stacking to the variance in the 
high-latitude sky. Last, we present an analytical approximation 
that relates the HRO between the filaments and J5 PO s t0 the E-B 
asymmetry. 

6.1. Stacked images of filtered maps 

Here we use the Planck D 353 , U 353 , £ 353 , and B 353 maps to 

analyse the £ range, between 30 and 300, over which we identify 
the filaments in the dust intensity map. Our choice of £ range 
has a large overlap with the angular scales, 40 < £ < 600, 
where the E-B asymmetry has been measured with the power 
spectra of dust polarisation. We stack all the bandpass-filtered 
maps, rotated by the mean orientation angle of the filament, 
as described in Sect. 5.1. The Q 353 and U ' 353 maps are the fil¬ 
tered Stokes 2353 and C /353 maps computed with respect to the 
axis of the filament. The mean stacked images of the bandpass- 
filtered maps are presented in Fig. 10. The sidelobes that appear 
in the images of (£> 353 ), ( 6353 )^ an d (^ 353 ) of Fig. 10 on both 
sides of the filament centre are coming from the filtering. The 
1 cr errorbar on the ( 2353 X (^ 353 )’ (^ 353 )’ an d (^ 353 ) images is 
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Fig. 11. Mean images of the Planck (^ 53 ^ 353 ^ (^ 353 ^ 353 )’ ^ 353 ^ 353 )’ an d ^ 353 ^ 353 ) ma P s over the A- filaments. 
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0.2 /a Kcmb, as computed from the difference of two polariza¬ 
tion HM maps. The filaments appear as a negative feature in the 
(Q 353 ) image and a positive feature in that of (£” 353 ). The (U' 353 ) 
and (# 353 ) images are consistent with a mean value zero. 

Next, we stack the products of two quantities, i.e., 

^ 353 ^ 353 ^’ ^ 353 ^ 353 ^’ ^ 353 ^ 353 ^ ^ 353 ^ 353 ^’ aiK * ^ 353 ^ 353 ^ 

We follow the same methodology as discussed in Sect. 5.1 
with oriented stacking at the centre of the filament. To avoid a 
noise bias in the square quantities, (£ 353 £ 353 ) and (£ 353 £ 353 ), 
we compute the cross-product of the two HalfMission maps 
(HM1 and HM2). For other quantities, (£ 353 £ 353 ), (^ 353 ^ 353 )’ 
and (£ 353 £ 353 ), we use the full-mission maps. The images pro¬ 
duced by stacking the products are presented in Fig. 11. Note 
that the maps shown in Figs. 10 and 11 are not used for data 
analysis, but still we comment on these images. We computed 
the mean stacked images from each of the two polarization DS 
or HR maps and all the features presented in Fig. 11 appear for 
both independent subset of Planck data. In the idealized descrip¬ 
tion of the filaments in Zaldarriaga (2001), we will expect the av¬ 
erage filament to appear similarly in the EE , ££, and TE maps. 
The differences between these three images shows that the real¬ 
ity of the dust sky is more complex than the idealized model. 


6.2. Measured E-B asymmetry 

We measure the E-B asymmetry at the filtering scale over the 
high-latitude (HL) region. We compute the variances (V) of the 
E 353 and 5 353 maps using the relations 


| Nhl 

V EE (HL) = Z ^353,HM1^353,HM2 = (46.6 ± 1.1)//K^ MB , 

(34) 

V BB (HL) = ^353,HM1^353,HM2 = (^9.1 ± 1.0)^K^ MB , 

(35) 


from the measurement at the power spectra level (Planck 
Collaboration Int. XXX 2014), probably because of the multi¬ 
pole range over which the ratio of the £353 and £353 variances 
are computed. 

We compute the covariances of the three maps over the high- 
latitude sky using the relations 


V TE (HL) = — ^ £353 £353 

V TB (HL) = — ^ £353 £353 
V EB (HL) = JJ— ^ £ 353 £ 353 


(124.1 ± 1.4)^K| mb , (38) 
(3.0±1.5) Ai K2 MB , (39) 

(-0.2±0.3)//K| mb . (40) 


The above covariances divided by V EE are listed in Table 1 . The 
ratio V TE /V EE that we find is consistent with the measurements 
of Planck Collaboration Int. XXX (2014) at the power spectrum 
level. 


6.3. Contribution of filaments to the variance of the E and B 
maps 

In this section, we compute the variance of the £ 353 and £ 353 
maps over the sky pixels used to produce the stacked images in 
Fig. 10, and compare the values with those measured over the 
HL region. These pixels are within the 7° x 5° patches, with an 
orientation angle 6 _, centred on the filaments (Sect. 5.1). These 
pixels define the grey regions in Fig. 12. We label them as SP 
and the rest of the high-latitude sky as O. The stacking procedure 
includes the filaments along with their surrounding background 
emission and, hence, effectively increases the selected fraction of 
the high-latitude sky, fi , from 2.2 % (filament pixels as described 
in Sect. 3.2) to 28 %. 

We compute the variance from the SP pixels using the rela¬ 
tion given in Eq. (34), 


where A/rl is the total number of pixels in the HL region. The 
ratio of the filtered £353 and £353 variances is 


V BB (HL) 
V EE (HL) 


0.62 ± 0.03 . 


(36) 


V EE (SP) = (137.5 ± 1.4)juK2 mb . (41) 

The sky variance of the £ 353 map in the high-latitude sky can be 
written as the sum of contributions from SP (/1 = 0.28) and O 
(1 - /1 = 0.72) regions. It is given by 


The uncertainty on V EE is computed by repeating the calcu¬ 
lation of Eq. (34) using the different cross-products, i.e., the 
two HalfRing (HR1 and HR2) and the two DetSet (DS1 and 
DS2) maps. We use the cross-HalfMissions as a reference for 
mean V EE . The 1 cr uncertainty on V EE comes from the differ¬ 
ences of these cross-products (DetSets minus HalfMissions and 
HalfRings minus HalfMissions). This 1 cr uncertainty is dom¬ 
inated by data systematics rather than statistical noise. This is 
in agreement with the uncertainties on power spectra over the 
same € range, as shown for / S k y = 0.5 in Figure 2 of Planck 
Collaboration Int. XXX (2014). The statistical noise on V EE 
is estimated from the cross-HalfMissions between their two 
HalfRing half-differences, N EE = (£353 HR1 - £ 353 ,^ 2)™!/ 2 x 

(^353,HRl ~~ ^353,HR2) HM2 / 2 aS 

CTyEE (HL) = = 0.02juK2 mb . (37) 

VMtL 

A similar procedure is applied to compute the uncertainty on 
V BB . The ratio of the £353 and £353 variances differs slightly 


V EE (HL) = /1 X V EE (SP) + (1 -/|)X V EE (O). (42) 


The ratio (£sp) of the variance from the stacked pixels to the 
total sky variance is given by 


D _ fix V EE (SP) 
SP V EE (HL) 


0.83 . 


(43) 


The value of £sp is expected to be high, since the filaments are 
bright structures on the sky. The pixels we used for stacking con¬ 
tribute 83 % of the total sky variance in the high-latitude sky. A 
similar result has been reported for the synchrotron emission, 
where bright filaments/shells are also measured to contribute 
most of the sky variance in polarization (Vidal et al. 2015). 

It has been noted that the structure in the Planck 353 GHz 
dust polarization maps is not fully accounted for by the filaments 
seen in the total dust intensity map. In particular, the local dis¬ 
persion of the polarization angle shows structures in the polar¬ 
ization maps that have no counterpart in total intensity (Planck 
Collaboration Int. XIX 2015). These structures are thought to 
trace morphology of £pos uncorrelated with matter structures 
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Fig. 12. Map of the selected pixels (grey colour) used in the 
stacking analysis. It covers 28 % of the high-latitude sky. Each 
tile in the image is a 7° x 5° patch around the filament centre and 
rotated by 6-. 


(Planck Collaboration Int. XX 2015). However, as our 5 sp value 
shows, these polarization structures do not contribute much to 
the variance of the dust polarization. 

In the same way as in Sect. 6.2, we compute the variances 
V EE , V BB , V TE , V TB , and V EB over the SP and O regions. Table 1 
presents the ratios of the variances computed over different sky 
regions. 


Table 1. The ratios of the variances computed from the selected 
pixels (SP) used in the stacking analysis and the rest (O) of the 
high-latitude sky (HL). 


Ratio 

SP 

O 

HL 

v BB /v EE 

0.66 ± 0.01 

0.51 + 0.05 

0.62 + 0.03 

yTE jy EE 

2.74 + 0.04 

2.48 + 0.15 

2.67 + 0.07 

yTB /yEE 

-0.07 + 0.04 

0.12 + 0.04 

0.06 + 0.03 

yEB/yEE 

-0.01 ± 0.02 

0.02 + 0.02 

0.00 + 0.01 


6.4. Analytical approximation 

Following Zaldarriaga (2001) and the description of the fila¬ 
ments in Eqs. (19) and (20), we can express the E- and 5-modes 
of a given filament as 

^353 ~ ~~ 2353 = ^353 C0S » (44) 

^353 ~ ““ ^353 = ““ ^353 S ^ n ’ (^ 5 ) 


where P 353 is the mean polarization intensity of the filament. 
The direct relation between E 353 and Q 35V and £35 3 and U' 353 
only holds for an idealised filament. For N idealized filaments 
oriented arbitrarily on the sky with respect to the GN, the ratio 
between the variances of the 5 and E maps is given by 


Y^_ <^ 53 * 353 > 

v EE ' <^53^53) " 2 ^ Cf w 2 ’ 


(46) 


which is expanded in terms of power spectra C BB and C EE un¬ 
der the assumption of statistical isotropy and homogeneity. The 
beam window function (we) is the filter function. Both the ob¬ 
served C BB and C EE dust power spectra follow a power-law 


model with the same slope a (Planck Collaboration Int. XXX 
2014). This reduces Eq. (46) to 

V BB £(2£+ l)A BB r a w] a bb 
V m - YSW + 1 ) ~F a vv) ~ A^E ■ (47) 

From the histogram of the upper panel of Fig. 6 , the distribu¬ 
tion of the angle A? - is known for our filament sample. Similar 
to the assumption made in Sect. 5.2, we assume that all the fila¬ 
ments have the same polarized intensity and therefore 


a bb (sin 2 2A| ) 

TFF *---= 0.66 

A EE (cos 2 2A| ) 


(48) 


We have computed the ratio A BB /A EE using the two independent 
subsets of the Planck data (HM maps) and find the same mean 
value of 0.66. This value of the A BB /A EE ratio based on this an¬ 
alytical model matches the observed mean value of 0.62 ± 0.03 
(Sect. 6.2). We note that the model value is directly inferred from 
the distribution of A|_- for our filament sample. If the HRO of 

A|_- was flat with uniform probability between -90° and +90°, 
we would have found equal variances in both the £353 and £353 
maps. 

In summary, we propose that the alignment between 5 P os 
and the filament orientations accounts for the E-B asymmetry 
in the range of angular scales 30 < £ < 300. The mean value 
of V BB /V EE for O region is consistent with the HE value within 
1.9 cr (Table 1). This shows that a similar alignment between the 
matter structures and 5 P os can be inferred over the rest of the 
high-latitude sky. Some high Galactic latitude sky areas, such as 
the BICEP2 field (BICEP2 Collaboration 2014), do not include 
any of the strong filaments from our study. The Planck 353 GHz 
polarization maps do not have the required signal-to-noise ratio 
to measure the A BB /A EE ratio for individual BICEP2-like fields 
(Planck Collaboration Int. XXX 2014). Therefore more sensitive 
observations will be needed to test whether our interpretation is 
relevant there. 


7. Relation to Galactic astrophysics 

In this section we place the paper results in the context of earlier 
studies about the filamentary structure of interstellar matter and 
its correlation with the Galactic magnetic field. 

Over the last decades, observations of interstellar gas and 
dust have been revealing the filamentary structure of the in¬ 
terstellar medium in increasing details. Before Planck and 
Herschel , the discovery of the infrared cirrus with the IRAS and 
Hi all-sky surveys was a main milestone in our perception of 
the structure of the diffuse ISM (Boulanger 1994; Kalberla & 
Kerp 2009). At high Galactic latitude, the Planck dust emission 
is tightly correlated with Hi emission at local velocities (Planck 
Collaboration Int. XVII 2014). In particular, all of the filaments 
in our sample have an Hi counterpart. They are selected on the 
353 GHz map but are seen at far-infrared wavelengths, in partic¬ 
ular the IRAS 100 pm map. 

The interstellar filaments seen in the Planck 353 GHz dust 
intensity map are not all straight. With our filament-finding al¬ 
gorithm, we have identified the straight segments with lengths 
L > 2°. Some of these segments are pieces of longer non-straight 
filaments. The 259 filaments in our sample make most of the sky 
variance in polarization as measured in our analysis. This is not 
a complete sample but other filaments at high Galactic latitude 
do not contribute much to the dust power in £-modes. 
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A number of studies, starting with the pioneering work of 
Goodman et al. (1990), have used the polarization of background 
starlight to investigate the relative orientation between the inter¬ 
stellar filaments and 2 ?pos- While most studies have targeted fil¬ 
aments identified in extinction maps of molecular clouds in the 
solar neighbourhood, a few studies have focussed on filaments 
seen in Hi emission (McClure-Griffiths et al. 2006; Clark et al. 
2014). These last two papers report a preferred alignment be¬ 
tween the filaments and 5 P qs i n th e diffuse ISM. The analysis of 
Planck data has complemented earlier studies providing greater 
statistics and sensitivity. Planck Collaboration Int. XXXII (2014) 
compare the orientations of matter structures identified in the 
Planck 353 GHz map with that of the Galactic magnetic field 
at intermediate latitudes. The alignment between the filaments 
and Z?pos reported in this paper becomes weaker for increas¬ 
ing column density (see Figure 15 of Planck Collaboration Int. 
XXXII 2014). Towards molecular clouds the relative orientation 
is observed to change progressively from preferentially parallel 
in areas with the lowest column density to preferentially per¬ 
pendicular in the areas with the highest column density (Planck 
Collaboration Int. XXXV 2015). The transition occurs at a col¬ 
umn density of 10 2L7 cm -2 . All the filaments considered in our 
analysis are much below this transition limit and observed, as ex¬ 
pected from earlier studies, to be statistically aligned with Z? PO s- 

Planck Collaboration Int. XXXII (2014) and Planck 
Collaboration Int. XXXV (2015) discuss these observational re¬ 
sults in light of MHD simulations, which quantify the respective 
roles of the magnetic field, turbulence, and gas self-gravity in the 
formation of structures in the magnetized ISM. The alignment 
between the magnetic field and matter structures in the diffuse 
ISM is thought to be a signature of turbulence. Simulations show 
that turbulent flows will tend to stretch gas condensations into 
sheets and filaments, which appear elongated in column den¬ 
sity maps (Hennebelle 2013). These structures will tend to be 
aligned with the magnetic field where the gas velocity is dynam¬ 
ically aligned with the field (Brandenburg & Lazarian 2013). 
Alignment also results from the fact that matter and the mag¬ 
netic field are stretched in the same direction because the field is 
frozen into matter. The change in relative orientation observed 
within molecular clouds might be a signature of the formation 
of gravitationally bound structures in the presence of a dynami¬ 
cally important magnetic field. Indeed, Soler et al. (2013) report 
a change in the relative orientation between matter structures and 
the magnetic field, from parallel to perpendicular, for gravita¬ 
tionally bound structures in MHD simulations. This change is 
most significant for their simulation with the highest magnetiza¬ 
tion. 

8. Conclusion 

We present a statistical study of the filamentary structure of the 
353 GHz Planck Stokes maps at high Galactic latitude, relevant 
to the study of dust emission as a polarization foreground to the 
CMB. The main results of our work are summarized as follows. 

We filter the intensity and polarization maps to isolate struc¬ 
tures over the range of angular scales where the E-B power 
asymmetry is observed. From a Hessian analysis of the Planck 
total dust intensity map at 353 GHz, we identify a sample of 259 
filaments in the high-latitude sky with lengths L > 2°. We mea¬ 
sure the mean orientation angle of each filament in this sample 
and find that the filaments are statistically aligned with the plane 
of the sky component of the magnetic field, B PO s, inferred from 
the polarization angles measured by Planck. We also find that the 
orientation of 5 P os is correlated with that of £ m ,pos in the solar 


neighbourhood. Our results show that the correlation between 
the structures of interstellar matter and 5 P os in the diffuse ISM 
reported in Planck Collaboration Int. XXXII (2014) for interme¬ 
diate Galactic latitudes also applies to the lower column density 
filaments (a few 10 19 cm -2 ) observed at high Galactic latitude. 

We present mean images of our filament sample in dust in¬ 
tensity and Stokes g 353 and C /353 with respect to the filament 
orientation (g ' 53 and U' 353 ), computed by stacking individual 
7° x 5° patches centred on each filament. The stacked images 
show that the contribution of the filaments is a negative feature 
with respect to the background in the Q' 353 image and is not seen 
in the U 353 image. This result directly follows from the fact that 
the histogram of relative orientation between the filaments and 
5 P os peaks and is symmetric around 0°. Combining the stacked 
images and the histogram, we estimate the mean polarization 
fraction of the filaments to be 11 %. 

We relate the E-B asymmetry discovered in the power spec¬ 
trum analysis of Planck 353 GHz polarization maps (Planck 
Collaboration Int. XXX 2014) to the alignment between the fil¬ 
aments and Z?pos in the diffuse ISM. The set of 7° x 5° patches 
we stack represents 28 % of the sky area at high Galactic lati¬ 
tude. The power of the 5-mode dust polarization computed over 
this area amounts to 83 % of the total dust polarization power 
in the high-latitude sky. We show with an analytical approxima¬ 
tion of the filaments (based on the work of Zaldarriaga 2001 ), 
that the HRO between the filaments and 5 P os may account for 
the Cf/Cf ratio measured over the high-latitude sky. Our in¬ 
terpretation could also apply to the E-B asymmetry reported for 
the synchrotron emission (Planck Collaboration X 2015), since 
there is also a correlation between the orientation angle of 5 P os 
and the filamentary structures of the synchrotron intensity map 
(Vidal et al. 2015; Planck Collaboration XXV 2015). 

Present models of the dust polarization sky (e.g., O’Dea 
et al. 2012; Delabrouille et al. 2013) produce an equal amount 
of power in E- and 5-modes for masks excluding the Galactic 
plane, because they ignore the correlation between the structure 
of the magnetic field and that of matter. Our work should mo¬ 
tivate a quantitative modelling of the polarized sky, which will 
take into account the observed correlations between the Galactic 
magnetic field and the structure of interstellar matter. 
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Appendix A: Hessian analysis 


In this appendix, we detail the implementation of the Hessian 
analysis to determine the orientations of the filaments. We start 
with the Z) 353 map at HEALPix resolution = 512. For each 
pixel on the sky, we estimate the first and second derivatives of 
D b 53 with respect to the Galactic longitude / and latitude b (as 
described in Monteserfn et al. 2005). The Hessian matrix of D 353 
is defined as 


H = 


#xx 

Hxy Hyy 


(A.l) 


where 




d 2 D b 

U ^353 

d 2 b 


d 2 D b 

tt _ 353 

xy cos b db dl ’ 
_ d 2 D b 53 
H yy cos 2 b d 2 l 


(A.2) 

(A.3) 

(A.4) 


We decompose the D b 53 map into a^ m coefficients using the 
“ianafast” routine of HEALPix and then use the “isynfast” 
routine in HEALPix to compute the second partial derivatives at 
each pixel. This method is computationally faster than the one 
used in Planck Collaboration Int. XXXII (2014). The smallest 
eigenvalue, /L, of the Hessian matrix is calculated as 


T_ 


2 (^x + Hyy 


-a). 


(A.5) 
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where a = ^l(H xx - H yy ) 2 + AH 2 y . The map of T_, displayed in 

the lower left panel of Fig. 2, highlights the filaments in the £> 353 
map. The filament orientation angle, 6 _, is calculated using the 
relation 


6 - = atan 




H yy “l - ct 


2 H, 


xy 


(A. 6 ) 


This follows the IAU convention, since it is measured from the 
GN and positive to the East direction. The formula for the fila¬ 
ment orientation angle 6 _ is equivalent to equation (9) of Planck 
Collaboration Int. XXXII (2014). 


S b 




A- [K C mb rad 2 ] 


Fig.A.l. Filtered realization of the Gaussian dust sky (upper 
panel) and the corresponding eigenvalue A- (middle panel). The 
distribution of the eigenvalue A- (lower panel) over the un¬ 
masked pixels is Gaussian. 


Appendix B: Gaussian realization of the dust sky 

Here we apply the SMAFF algorithm on a Gaussian realization 
of the dust total intensity map using the same methodology as 
described in Sect. 3.2. The power spectrum is modelled as C BT oc 
l a , where a is the slope of the power law. We choose a = -2.4 
in the multipole range from l - 30 to 300 (Planck Collaboration 
Int. XXII 2015). 

We compute a Gaussian realization of a sky map on a 
HEALPix grid with an t~ 1A power spectrum, and use the spline 
wavelet transform to filter it. Hereafter, we refer to the filtered 
simulated Gaussian map as the “ S b map”. We compute the 
smallest eigenvalue T_ of the S h map using the Hessian analysis 
described in Appendix A. The S h map and its corresponding A- 
map are shown in the upper and middle panels of Fig. A.l. 

We only consider the high-latitude sky. With the upper 
threshold K- = m A _ - 3 cr A , we remove pixels as shown in the 
lower panel of Fig. A.l. Running our friend-of-friend algorithm 
on the S b map, we do not detect any filament with L > 2°. For 
higher values of the threshold X_, we detect such filaments in 
the S h map. We call them “weak” filaments. The threshold fac¬ 
tor X_, used in our study, is a key factor to separate the strong 
filaments from the weak ones. The main results of this paper fol¬ 
low from the statistical properties of strong filaments. However, 
in Appendix C we demonstrate that they still hold for the weak 
filaments. 


Appendix C: Effect of A- thresholding on the 
filament count 

We apply the threshold = m A _ - 3 cr A in Sect. 3.2 to find 
strong filaments in the D b 53 map. In this section, we change 
the threshold K_ to quantify its effect on the HRO between the 
filaments and 2 ?pos- We choose different thresholds and divide 
the selected A- coherent structures into two categories, namely 
strong and weak filaments. The strong filaments are selected 
from the sky pixels that have A- < m A - 3 cr A _. Whereas, 
the weak filaments are selected from the sky pixels that have 
m A _ - 3 cr A _ < A- < m A _ - p cr A _ , where p is a factor to defined 
the threshold = m A _ - p cr A _ . In Table C. 1 , we list the number 
of weak filaments for different values of p. 

By choosing a threshold = m A _ - 1 cr A _ and running our 
friend-of-friend algorithm on the D b 53 map, we double the num¬ 
ber of filaments with L > 2°. For this larger set, we compute the 
HRO between the filaments and Bpos, and present it in Fig. C.l. 
The 1 cr dispersion of the angle difference remains the same as 
for our nominal set. This shows that the alignment between the 
structures of interstellar matter and Bp os holds even when in¬ 
cluding the weak filaments. 

Including the weak filaments and their surrounding back¬ 
ground emission in the analysis of Sect. 6.3 effectively increases 
the sky fraction from 28 % to 50 % of the high-latitude sky. With 
more sky coverage, the ratio of the variances, i.e, V BB /V EE , 

yTE/yEE ? yTB /V EE ? and yEB/yEE i§ dose to the HL region val _ 

ues quoted in Table 1 . 

Appendix D: Uncertainties on the angles 6- and A 

In Sect. 4 we present our statistical analysis of the relative ori¬ 
entation between the filaments and 2 ?pos • The filament orienta¬ 
tion angle, (L, is a single number computed from the D b 53 map. 
To estimate the error bar on 0_, we propagate the noise in the 
£>353 map (N D353 ) through the Hessian analysis. We simulate a 
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Table C.l. Total filament count, including the strong and weak 
filaments, as a function of the threshold K- on the A- map de¬ 
rived from the map. 


Factor 

P 

m A _ 

Filament count 

weak 

- 3 cr A _ < A_ < m A _ - p cr A _ /L 

strong 

<m A _ -3 cr A _ 

total 

3 

0 

259 

259 

2 

53 

259 

312 

1.5 

140 

259 

399 

1 

259 

259 

518 



Fig. C.l. Same as the upper panel of Fig. 6 , when including an¬ 
gle differences from the weak filaments. The width of the 
histogram does not change. 


Gaussian realization of the N Ds53 map and add it to the £>353 map. 
We filter the noise added to the £>353 map in the same manner as 
the £>353 map. For each filament in our sample, we compute the 
change of its orientation angle from the added noise using the 
formula 

A0_ = 6 _ (D b m + N b h J - 0_ (D% 3 ). (D.l) 

The distribution of the angle difference is presented with red 
colour in Fig. D. 1. The mean of the distribution is 0° and the 1 cr 
dispersion is 0.2°. This shows that the noise in the £>353 map has 
little impact on our estimate of the filament orientation angle. 

Next, we compute the uncertainty on the orientation angle of 
2 ?pos- We use the two HM Planck polarization maps, and com¬ 
pute the orientation angle of #pos using Eq. (2). For each fila¬ 
ment in our sample, we compute the difference between the two 
values obtained from the two HM maps, 

AT = 0.5 x [T H mi ~^hm2]- (D-2) 

The histogram of the angle difference is presented with blue 
colour in Fig. D.l. The mean of the distribution is 0° and the 
1 cr dispersion is 2°. The error on the orientation angle of #pos 
is small compared to the width of the HROs in Fig. 6 . 

Last, we assess the impact of the map filtering on our anal¬ 
ysis, computing the mean orientations of the filaments, over the 
same set of selected pixels, on the £>353 map. We compute the 


angle difference of the filament orientation derived from the fil¬ 
tered and unfiltered Planck 353 GHz maps 

A© = 9- (D 353 ) - 9. (D b 353 ). (D.3) 

The histogram of A0 is plotted with green colour in Fig. D.l. 
The mean of the distribution is 0° and the 1 cr dispersion is 7.0°. 
All the uncertainties measured in 6_ and X are small compared 
to the dispersion of the HRO measured in Fig. 6 . This means that 
the data noise (for the intensity and polarization maps) and the 
filtering of the data are not critical for our study. 



°-30 -24 -18 -12 -6 0 6 12 18 24 30 

[deg] 

Fig. D.l. Histogram of changes in (red) and A (blue) due to 
noise present in the Planck data. The histogram of changes in 
6_ from the difference of the filtered and unfiltered £>353 maps is 
shown in green. 
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